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ABSTRACT 

Since few decades, asteroseismology, the study of stellar oscillations, enables us to probe the interiors 
of stars with great precision. It allows stringent tests of stellar models and can provide accurate 
radii, masses and ages for individual stars. Of particular interest are the mixed modes that occur in 
subgiant solar-like stars since they can place very strong constraints on stellar ages. Here we measure 
the characteristics of the mixed modes, particularly the coupling strength, using a grid of stellar 
models for stars with masses between 0.9 and 1.5 Mq. We show that the coupling strength of the 
£ = 1 mixed modes is predominantly a function of stellar mass and appears to be independent of 
metallicity. This should allow an accurate mass evaluation, further increasing the usefulness of mixed 
modes in subgiants as asteroseismic tools. 

Subject headings: stars: oscillations, stars: interiors, methods: data analysis 



1. INTRODUCTION 

Asteroseismology allows stringent tests of stellar mod- 
els and can provide accurate fundamental properties 
of individual stars. For solar-type stars on the main 
sequence, the observed oscillations are p modes, for 
which the restoring force arises from the pressure gra- 
dient. These are approximately regularly spaced in fre- 
quency, fohow^ngclosjsly the so-called asymptotic re - 
lation (jVandakurovi 119671 : iTassoull 119801: iGoughl 119861 ). 
However, the oscillations of post-main-sequence stars 
show departures from this regularity that are due to the 
presence of mixed modes. 

Mixed modes have p-mode character in the stellar en- 
velope and g-mode character in the core. They occur 
in evolved stars (subgiants and red giants), in which the 
large density gradient outside the core effectively divides 
the star into two coupled cavities. This leads to mode 
bumping, in which mode frequencies are shifted from 
their regular spacing and no longer follow the asymptotic 
relation. Mode bumping in s ubgiant stars was first ob- 
served and modeled in r\ Boo (Kields en et al.111 995, 2003; 
IChristensen-Dalsgaard et al.l 119951: ICarrier et al.l I2005D 
and 3 Hvi (jBedding et all 120071: Brandao et all 120111 ). 
More recently, asteroseismic space missions have pro- 
duced man y more examples, inclu ding the CoRoT target 
HD 49835 (IDeheuvels etHllMol) an d a growing number 



Metcalfe et aLll2010t iMathur et al.1 



of Kepler stars (e.g , . 
IMllCampante et all 1 2 6T1 

Mixed modes carry valuable information on the in- 
ternal structure and evolutionary state of stars. Their 
frequencies change quickly wit h time as they underg o 
avoided crossings (jOsakil il975t lAizenman etaTJ 119771) . 
potentially providing stellar ages with a precision down 
to a few Myr, or a relative uncertainty of « 1% 
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HMetcalfe et al.ll2010h . 

Mixed modes arise from a resonant coupling between p 
and g modes and can be well represented by a system 
of coupled oscillators. In this Letter, we use such a rep- 
resentation to model their frequencies in subgiant stars 
and determine the coupling strength between the modes. 
We show that the coupling strength depends strongly on 
stellar mass but only weakly (or not at all) on metallicity, 
hence lifting the often problematic degeneracy between 
those two variables. 

2. A MODEL FOR MIXED MODES 

Mixed modes in stars occur when the p and g modes 
are couple d, which leads to avo i ded cro ssings and mode 
bumping. IDeheuvels fc Michel (|2010al) suggested that 
an avoided crossing in a subgiant star can be well- 
represented by a system of (n — 1) p-mode oscillators, 
each coupled with a single g mode. This can be modeled 
by a system of n differential equations, one for each of 
the p modes and one for the g mode: 



d 2 yx(t ) 
dt 2 



(1) 



d 2 y n - 1 (t) 
dt 2 

d 2 y n {t) 
dt 2 



= -U^y n + OLX,nV\ + ■•• + an-l,nl/n-l 



The yi(t) terms are the displacements of the modes and 
ai.n are the coupling coefficients between the ith p mode 
and the g mode. The frequencies ui 7Ti and w 7 correspond 
to fictitious, pure p and g modes that would exist if their 
cavities were not coupled. In or der to avoid ambiguity, 
we follow lAizenman et all (|1977l ) by r eferring to these as 
7r and 7 modes, respectively (see also Bedding 2011a). 

We seek oscillatory solutions of Eq.[T]with angular fre- 
quencies VL — {oji, UJ2, ...,w„}, which means we need to 
solve the following system: 



AY = tt z Y, 



(2) 
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where Y = {y 1 ,y 2 , ...,y n } and 
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(3) 



Following [Deheuvels k Michell (|2010aD . we assumed the 
coupling a,i ;n = a to be the same between all modes. 
This is a reasonable assumption provided a varies slowly 
with frequency, which turns out to be the case in the 
models we have studied. 

In more evolved stars (red giants), the situation is re- 
versed and several g-modes may be coupled to a sin- 
gle jP - mode dDupret et al.l l2009t | Christensen-Dalsgaardl 
120111: lBeddmdl2011at iStellol 120111) . Between these two 
extremes, one can consider inter mediate cases with mu l- 
tiple g-modes and p-modes (e.g.. Idi Mauro et al.ll2011[ ). 
whose power spectra may become hard to interpret. 

For the case of n = 2, which is a single 7r mode coupling 
to a 7 mode, Eq. [2]can be solved analytical ly to give a 
(Do 



pair of solutions (|Deheuvels k Michell[20T0ah : 



± 



4a 2 . 



(4) 



The black curves in Figured] s hows these two solutions 
in a replicated echelle diagram (|Beddindl20rTbl) . 

We note that the coefficient a is not an ideal measure 
of the coupling strength between 7r and 7 modes because 
it varies with frequency and is therefore different for each 
avoided crossing. A better measure is the frequency sep- 
aration of the two solutions at the avoided crossing. That 
is, we focus on the point of resonance by setting w 7 = uj v 
and measure the frequency difference between the two 
solutions: 



Sujj = uj+ — uj- 



It follows from Eq.|4] that 



or, in cyclic frequency, 



SiSy = 5ujy/2ir 



(5) 



(6) 



(7) 



For the case of n > 2, which involves several tt modes 
coupling to a 7 mode, Eq.[2] can be solved numerically. 
We find a set of solutions, shown by the red curves in 
FigureQ]for the case of n = 8. The value of 6u 7 is shown 
by the solid red line. 

The minimal separation is almost independent of the 
number of interacting modes (black and red lines are 
close to each other at ). The stronger the cou- 

pling, the larger the separation. 

3. FITTING METHOD 

For each avoided crossing we have a sequence of ob- 
served frequencies that correspond to the mixed modes. 
Each observed frequency must satisfy Eq.[2] Provided we 
have enough observed frequencies, we can determine the 
elements of the matrix A, which gives the tt and 7 mode 
frequencies, and the coupling strength a. 

To carry out this fitting process robustly, we ap- 
ply a Bayesian approach (Maximum a Posteriori Ap- 
proach, hereafter MAP). Assuming a likelihood function 




Figure 1. Replicated echelle diagram showing the effect of the 
coupling between the tt and 7 modes in two configurations: two 
coupled modes (black), and several coupled modes (red) with tt- 
modes frequencies spaced by Au. u~/ is arbitrarily fixed to 10. 
Where the tt and 7 modes cross, the frequency separation is <5i/ 7 
(diagonal red line). Su 7 is almost independent of n, the number of 
interacting modes. 



L(^obs|^4.), one can regularize this function by penalizing 
it with a quantity p(A), called the prior, where f b s rep- 
resents the observed frequencies (the data). Using Bayes' 
theorem, we write the statistical criteria defining the fit 
quality that we seek to maximize as 

ln[p(AK bs )] = ln[L(u ohs \A)} + ln[p(A)] - ln[C]. (8) 

Here, C is a normalization constant and p{A\u Q i os ) is the 
posterior probability (the probability of A given v hs)- 
The linear system in Eq.[2] must be solved at each it- 
eration of the maximization process in order to com- 
pare the observed frequencies ^ b s with the calculated 
frequencies v ca \ c - The expression we choose to minimize 
is \ J2i C^obsCn) - i^caic (n)) 2 /a- 2 . Thus the log-likelihood 
is, 



ln[L(i/ obs |A)] cx 



^caic(n)) 2 /cr 2 . (9) 



where a n — 1 for model frequencies. 

Having defined the likelihood, we must choose expres- 
sions for the priors. No priors on a and were applied 
and our attention will focus on . The pseudo-modes 
7r and 7 are expected to behave as pure p and g modes. 
To a good approxim ation, the 7r-m odes should follow the 
asymptotic relation ([Tassoul 198(3): 



v w (n) 



(n+ - +e)Ai>- 



(10) 



where n is the radial order, e is an offset associated with 
stellar surface effects, is the large separation (related 
to the mean density of the star), and 5vqi is sensitive to 
the core properties. In this Letter, we only consider the 
case of I = 1 (dipole modes). 

Discontinuities inside stars introduce frequency oscil- 
lations as a function of n, not included in the asymptotic 
relation. These oscillations are not expected to be sig- 
nificant, say less than a few percent of Av. Thus the 
7r-mode frequency variations are expected to be smooth, 
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as seen in the echelle diagram for the p-modes of main se- 
quence stars. A way to satisfy this condition is to impose 
a smoothness condition on some pth derivative terms of 
Vn(n). We note that the second derivative in n of Eq.fTOl 

denoted A 2 u{n), leads to A 2 v{n) = w 0. Thus one 
may limit local strong deviations from a regular pattern 
by imposing a Gaussian prior on A 2 v(n), 



p(A 2 v{n)) 



1 



27T<TA 2 i/ 



exp 



A 2 v(n) 



(11) 



Here, a a?v plays the role of a relaxation constraint and 
must be chosen to ensure enough freedom, but not too 
much, in order to efficiently smooth the frequency pro- 
file. A trial-and-error procedure showed that o ap-v ss 2 
fiHz offers a good compromise. With such a smoothness 
condition, the 7r-modc deviation from a strictly regular 
pattern of frequencies is locally described by a second- 
order polynomial function of n, and the solution belongs 
to the family of spline functions. The smoothness con- 
dition acts locally and does not restrict the 7r-modes to 
follow the asymptotic relation, globally. 

We define an additional (global) condition, in order 
to avoid strong departures from Eq|TU] We expect the 
£ = modes and I = 1 it modes to be distributed along 
two parallel ridges in the echelle diagram. Hence, the 
large separation Av of the £ = mode, is expected to 
be approximately equal to the large separation Av n of 
the £ = 1 7r modes. Thus a second Gaussian prior was 
applied to the quantity Av n — Av, with ctav^ ~ l/y/N 
/xHz, TV being the number of modes. 

4. RESULTS FOR STELLAR MODELS 

We have applied our method to frequencies of £ = 1 
modes calculated from stellar models. Mo dels were se- 
lected from the grid described in detail by iStello et"aL1 



2009) . This grid was genera ted with the ASTEC code 
Christcnscn-Dalsgaard 2008 1) using the simpl e but fast 



EFF equation of state (jEggleton et al.1 119731 ) , a fixed 
mixed length parameter at q;mlt = 1.8, and an initial 
hydrogen abundance of X = 0.7 . The opacities were cal- 
culated using the solar mi xture dGrev cssc & Noels| [T993() 
and the opacit y tables of iRogers fc Iglesiasl (|1995|) and 
iKurucd ([1991) . Rotation, overshooting and diffusion 
were not included. 

Models with masses in the range 0.9 to 1.5 M were 
explored with a step size of 0.1 Mq. All models were 
computed for solar metallicity (Z — 0.017). In addition, 
lower and higher metallicities (Z = 0.011 and Z = 0.028) 
were considered for a subset of masses (1.0 Mq, 1.2 Mq, 
1.3 M q and 1.5 Mq). 

In general, stars have mixed 1=1 modes while still on 
the main sequence. However, these occur at low frequen- 
cies that lie outside the envelope of observable modes. 
Mode bumping only starts to be detectable once a star 
has entered the subgiant phase. As discussed in the In- 
troduction, each avoided crossing is associated with a 
g m ode trapped in the core of t he star (a 7 mo de; see 
also IDeheuvels fc Michell I2010M IBeddind [20ilal) . The 
7-mode frequencies increase as the star evolves and we 
tracked each one. As additional avoided crossings ap- 
peared at low frequency, they were also incorporated into 
Eq.[2j In this way, the stellar model was followed dur- 



ing a significant part of the subgiant phase, as shown in 
the HR diagram in Figure[2j Towards the end of this 
phase, the frequency patterns became too complex (too 
few modes between each avoided crossing) and a stable 
fit was difficult to obtain. We secured reliable results for 
models with up to three observable avoided crossings, 
each of them characterized by its frequency v~, and its 
coupling strength (measured by bv^ — see Eq.[7]). 

Figure [3] shows the results, with each panel showing the 
coupling strength as a function of the avoided crossing 
frequency. The central frequency of the modes v max de- 
creases over time and thus, goes from right to left along 
a star's evolution in these diagrams. The systematic er- 
ror on the determination of bv^ is about 1.5 /iHz. This 
uncertainty arises from Eq.[2] not being a perfect repre- 
sentation of the behavior of the frequencies of the stellar 
models. In Figured we show only the highest-frequency 
avoided crossing. The different tracks show models with 
four different masses and three different metallicities. In- 
terestingly, the co upling strength depends st rongly on 
mass (as noted by IDeheuvels fc Michel 1201 lh but only 
weakly on metallicity, if at all. The systematic error pre- 
vents us from assessing whether the small variations of 
the coupling strength are really due to the metallicity. it 
suggests that we may be able to use the coupling strength 
to break the degeneracy between mass and metallicity, 
which is often a problem in asteroseismology. 

As mentioned above, 7 modes frequencies increase as 
a star evolves and several avoided crossing may enter the 
range of observable frequencies. In Figured we show 5v~ 
for the first three avoided crossings for evolving models 
with four different masses. They all have solar metal- 
licity. For a given mass, all three curves closely follow 
the same path in the diagram (but with a time delay). 
This confirms that b IS-y , CIS defined by Eq.[?l is an excel- 
lent measure of the coupling strength between the p- and 
g-mode cavities. 

In Figure^ and [HJd we see rapid changes in the cou- 
pling strength at the low-frequency end. This reflects 
rapid changes in the extent of the evanescent zone that 
separates the inner g-mode and outer p-mode cavities. 
Indeed, the coupling strength depends on the extension 
of this zone, in the sense that a smaller evanescent zone 
leads to stronger coupling. In Figure^ the tracks for the 
different avoided crossings (for a given mass) no longer 
follow a single path at low frequency. During the tran- 
sition between the main sequence and the post main se- 
quence stage, the density of the star and thus, the shape 
of the evanescent zone, vary on very short timescale. 
Thus the evanescent zone has changed significantly be- 
tween the times when the first and third avoided crossing 
pass the same frequency. For the most evolved stars, 
which means when several 7 modes have a frequency 
close to v m ax, the diagram is also problematic, as the 
low-frequency region of the spectrum is densely popu- 
lated with 7 modes, making it hard to obtain a stable 
fit. 

To summarize, we conclude that the coupling strength 
of mixed £ = 1 modes, as measured by the minimal sep- 
aration bvy, is a useful observable in subgiant stars, pro- 
vided they are not too evolved. In particular, bv-y de- 
pends mainly on the stellar mass and is almost indepen- 
dent of metallicity and of which avoided crossing is being 
measured. We will now use this result to estimate the 
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Figure 2. HR diagram for models with M = 0.9 — 1.5Mq at solar 
metallicity. The thick colored lines show the range over which the 
three highest frequency avoided crossings were tracked. 

masses of subgiants for which observed frequencies have 
been published. 

5. APPLICATION TO OBSERVED FREQUENCIES 

We have applied the approach described above to the 
seven subgiants listed in Table [1] We use the published 
frequencies, which were measured either by ground-based 
spectroscopy, or from space by the CoRoT or Kepler mis- 
sions. 

Figure[5J: shows the results. The symbols show the re- 
sults of applying our method to the seven subgiants. The 
colored tracks are from the analysis of theoretical models 
discussed above, for a wide range of masses and a single 
metallicity, and only for the highest-frequency avoided 
crossing (the 71 mode). Crosses along each curve are 
spaced by 150 Myrs. We see that the 7-mode frequen- 
cies may vary by several hundreds of /iHz within this time 
scale while the precision of measure of these frequencies 
is only of few /iHz. Thus 7-mode frequencies provide a 
stringent constraint on stellar age. 

We note that any potential deviation between the 
tracks of different metallicity (Figure^) and for different 
avoided crossings (Figure[3jD) are within the uncertainties 
with which we can measure Svy. We can therefore super- 
impose results from more than one avoided crossing per 
star on the same diagram (Figured]:). The quality of the 
data allowed us to extract accurately two avoided cross- 
ing in three Kepler stars (Mulder, Boogie and Gemma). 

An example of a fit is shown in Figure|4] for Mulder. 
Despite its simplicity, the linear system of Eq.[5] can be 
useful to predict the position of previously unidentified 
mixed modes. Indeed, our method predicts an extra 
t = 1 mixed mode, hardly seen on the power spectrum 
becaus e it lies close to an I = 2 mode (see iBeddingl 
l2011bl for another example). This demonstrates that our 
method would allow us to identify I = 1 mixed modes in 
complicated power spectra. 

The results of the fit are given in Table [TJ The first 
three columns give the identification number, the name 
and the observed ^ max of each star. Columns 4-7 give 
the fitted values of za, and Su~ for the avoided crossings. 
The last two columns give our inferred stellar mass, and 
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Figure 3. (a): (5f 7 as a function of the 7-mode frequency for i/ 7 i 
and for I.OMq (red lines), 1.2Mq (green lines) and 1.3Mq (blue 
lines). Dashed, solid and dotted lines are for three metallicitics 
Z.(b): Svj as a function of v 1 i (solid lines), i/ 7 2 (dotted lines) and 
Vj3 (dashed lines) at Z = 0.017. (c): 8u-y as a function of u~/i 
for stellar masses M = 0.9 to A/ Q = 1.5, at Z = 0.017. Crosses 
are spaced by 150 Myr. Symbols correspond to the seven ana- 
lyzed stars. Whenever reliable, fit values for the second measured 
avoided crossing are shown (gray symbols). Doted lines indicate 
when v 7 i = 0.75^ ma x, f 7 i = fmax and v 7 i = 1.25t/ m ax- They 
indicate the best region to observe the first avoided crossing. 
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is the same for all avoided crossings and for a given star. 

In this Letter, the effects of varying the mixing length 
parameter and the initial helium abundance have not 
been tested. They are known to be correlated to stellar 
mass and deserve further investigation. 



Figure 4. Replicated echelle diagram (Bedding 2011b) for Mulder 
Kepler target). Two avoided crossings are visible. The greyscale 
shows the observed po wer spectrum and open symbols show the 
frequencies reported by Campantc ct al. (2011). The model pre- 
dicts an extra mixed mode (circles), not found by Campantc ct al. 
(2011) because it lies on an £ = 2 mode. 

the published mass found using conventional asteroseis- 
mic analysis. For the stars with published masses, our 
results agree well. However, in the case of HD 49835, 
the precision on the determination of bv~ j is low (and so 
for the mass) because the avoided crossing occurs right 
at the lower boundary of the observed modes, with only 
one i — 1 mode below i> 11 . 

The uncertainties on our mass estimates were ob- 
tained by linear interpolation between the tracks of dif- 
ferent masses (cf. F igureBl :). As already noticed by 
IDeheuvels fc Michell (poTTf T mass and coupling are in- 
versely proportional. Thus, uncertainties increase with 
mass because high-mass tracks are closer together than 
low-mass tracks. In the best cases, the total uncertainty 
(quadratic sum of the uncertainty of the measurement 
and of the systematic error) is a few percent, and in 
other cases it may be as large as 10%. 

The results presented here are based on observed I — 1 
mixed modes. However, the mean large separation com- 
puted from the t — frequencies provides an indepen- 
dent way to measure stellar masses. Thus, the global 
accuracy, can be increased by using inferred masses from 
both the coupling strength and the mean large separa- 
tion. 

With observational data, not all the mixed modes can 
be seen because the mode envelope has a limited width. 
In order to indicate this in Figure|3];, the upward-sloping 
curves show where v lx = 0.75z/ max , v lx — v max and 
v 11 = 1.25^ max . Thus, these curves show the most likely 
region in which to observe the first avoided crossing. 
The most likely region to observe the second and third 
avoided crossings can be obtained by shifting these lines 
towards lower frequencies by w 200/xHz and « 400/^Hz, 
respectively. 

Theoretically, the coupling strength should depend on 
the extent of the evanescent zone between the g-mode 
cavity in the core and the p-mode cavity of the envelope. 
Thus, it would be interesting to understand the exact 
nature of this relation, and why the minimal separation 
8v 1 seems to depend strongly on stellar mass whereas it 



6 



Benomar et al. 



Table 1 

Measured mixed modes properties (minimal separation <5i/ 7 and 7-mode frequency, v 7 ) and inferred mass M (followed by its uncertainty and 
systematic error). For comparison, Mi; t lists published masses derived by traditional seismic model fitting. 



Id. Number 




Nickname 


^max (jjBz) 


Uj 1 (f^Hz) 


8v 11 (/iHz) 


i/ 72 (/xHz) 


SUj 2 (/iHz) 


M (M ) 


M Ut (Mo) 


KIC 10273246 


(a) 


Mulder 


842 


970.7 ±3.2 


24.6 ±1.4 


707.3 ± 1.2 


31.4 ±0.9 


1.16 ±0.03 (±0.04) 




KIC 10920273 


(a) 


Scully 


974 


854.1 ± 2.6 


30.6 ±2.3 






1.04 ±0.04 (±0.04) 




KIC 11234888 


(b) 


Tigger 


675 


718.4 ±5.6 


23.5 ±1.3 






1.27 ±0.03 (±0.10) 


1.33 ±0.26 


KIC 11395018 


(b) 


Boogie 


830 


970.2 ± 10.6 


23.3 ±1.7 


747.7 ±0.6 


25.6 ±0.9 


1.21 ±0.06 (±0.04) 


1.25 ±0.24 


KIC 11026764 


(c) 


Gemma 


857 


920.7 ±3.1 


25.1 ± 1.0 


710.2 ±5.3 


26.3 ±2.1 


1.17 ±0.03 (±0.05) 


~ 1.13 or 1.23 


HD 49835 < d > 






1013 


758.3 ± 3.2 


22.1 ±3.2 






1.30 ±0.11 (±0.06) 


1.25 ±0.05 


P Hydri < c ) 






960 


975.5 ± 1.6 


29.3 ±1.2 






1.02 ± 0.06 (±0.04) 


1.08 ± 0.03 



M ICampante et al.l SWHh 
(b) IMathur et all d20T!fi 
< c ) IMetcalfe et al.l |[201CT ) 

( d ) IDeheuvels eTall poToT ) and IDeheuvels fc Micheil pOlT l 

( e ) IBrandao et"all ll20Tll ') 
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